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Experimental and numerical studies were carried out to evaluate the effect of vortex generator on a small cylin- 


drical protrusion at Mach number 2.0. The experiments were performed using the supersonic blow down wind 


tunnel on different heights of cylindrical protrusion with vortex generator placed ahead of them. The upstream 


and downstream flow around the cylindrical protrusion is influenced by vortex generator as is observed using 


both visualization and pressure measurement techniques. Numerical studies using three dimensional steady im- 


plicit formulations with standard k-œ turbulence model was performed. Results obtained through the present 


computation are compared with the experimental results at Mach 2.0. Good agreements between computation and 


experimental results have been achieved. The results indicate that the aerodynamic drag acting on cylindrical pro- 


trusion can be reduced by adopting vortex generator. 


Keywords: Cylindrical protrusion, Vortex generator, Shock wave boundary layer interaction, Drag control 


Introduction 


High speed flow around any protruded obstacles leads 
to a Shock Boundary Layer Interaction (SWBLI). Even 
though the design version is streamlined, the production 
versions of any aerospace vehicle always consist of sev- 
eral kinds of protrusions all over their external surface. 
The presence of these protrusions can alter the surface 
pressure distribution and generate adverse pressure gra- 
dient due to shock boundary layer interaction. Due to the 
presence of this adverse pressure gradient ahead of pro- 
trusion, the incoming boundary layer may separate lead- 
ing to a complex phenomenon in the vicinity of protru- 
sion. The bow shock contributes towards the drag of the 
whole body and the Shock Wave Boundary Layer Inte- 
raction can lead to larger structural loading on adjacent 
surfaces. Although many researchers contributed towards 
different geometries of protrusion, the height effect was 
not always considered as one of the main geometric pa- 
rameters. In case where the height of protrusion is very 
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small or is of the same order of local boundary layer thic- 
kness, the effect of height of the protrusion in the vicinity 
of flow field cannot be neglected. Westkaemper [1] pro- 
posed an empirical relation to obtain the separation dis- 
tance ahead of a cylindrical protrusion for different di- 
ameter and heights at different Mach numbers. Sedney [2, 3] 
performed extensive flow visualization studies and re- 
ported the results of different size of cylindrical protru- 
sions. 

The extent of upstream separation distance is mainly 
caused by the height of protrusion, the incoming boun- 
dary layer thickness (6) and the freestream Mach number 
[4, 5]. The effect of short protuberances on interactive 
flow field at hypersonic speed was studied by Su-xun Li 
et.al [6]. Interaction of three dimensional protuberances 
with a supersonic turbulent boundary layer was studied 
by extensive numerical simulation with Hybrid RANS/ 
LES model by A.Frendi [7]. To control this shock wave 
boundary layer interaction flow field which also occurs 
in supersonic inlets and normal shock interaction with 
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Nomenclature 
Ap vortex generator incidence angle c vortex generator chord length (mm) 
Cp coefficient of pressure y+ non-dimensional distance 
Cp coefficient of drag 
D diameter of cylindrical protrusion (mm) Greek letters 
H, height of cylindrical protrusion (mm) Ay first cell distance 
Hyg height of vortex generator (mm) 0) specific dissipation rate (m’/s*) 
Mə» Freestream Mach number ò boundary layer thickness (mm) 
x stream-wise direction measured from cylinder 

center Subscripts 
Y normal direction measured from cylinder center vg vort Gx generator 

c cylinder 

Z span-wise direction measured from cylinder center o0 freestream 


the boundary layer, several active and passive methods 
have been employed like bleeding of boundary layer, 
adopting vortex generator (VG) etc. VG is one of the 
physically robust passive methods which creates stream- 
wise counter rotating vortices to energize the boundary 
layer and mainly helps in delaying the separation. 

But there is also a possibility of reducing the aerody- 
namic drag by altering the flow field around the protru- 
sion. The smaller VGs are characterized by heights less 
than or equal to the flow boundary layer thickness and 
have been referred to as “low -profile” or “micro” VGs, 
particularly, micro-vanes and micro-ramps. Control of 
SWBLI using an array of micro-vanes and micro-ramps 
was reported by Anderson et al. [8]. He also demon- 
strates the Response surface Methodology (RSM) to de- 
termine the optimal designs of micro-array actuation for 
controlling the SWBLI. Experimental surface separation 
pattern is obtained by Babinsky [9], in which a series of 
distinct accumulated oil lines mark the location of sepa- 
ration lines. Using his experimental results of flow sepa- 
ration, a vortex model was derived. The vortex system 
around VG consists of 4 pairs of vortices, i.e., the horse- 
shoe vortex due to the leading edge separation of VG, the 
primary vortex by VG, two secondary vortices under the 
primary vortex - one lies on the bottom of the plate and 
the other lies at the side of VG. Lu et al. [10] reported 
experimental and numerical investigation on several 
types of micro-VG at supersonic speeds. 

Detailed study of flow interactions with micro-ramps 
on a supersonic boundary layer at M=3.0 was investi- 
gated using monotone integrated Large Eddy Simulation 
(MILES) and Reynolds Averaged Navier-Stokes (RANS) 
by Loth and Lee [11]. Ramp induced interaction studies 
were carried out using micro-vortex generators by Pierce 
et al. [12] who studied the flow features of micro-VG 
with l-inch diameter circular cylinder by conducting 
fluorescence surface flow visualization. Galbraith [13] 
investigated on the control of shock wave boundary layer 


interaction using multi row micro-ramp actuators and it 
was able to improve the boundary layer health down- 
stream of the separation. 

In the present study, an attempt has been made to util- 
ize the advantage of vortex generator to avoid flow sepa- 
ration due to shock wave boundary layer interaction 
ahead of the protrusion with a possibility of reducing the 
pressure drag due to protrusion. The effect of protrusion 
height on the flow features around a cylindrical protru- 
sion is computationally investigated. 


Experimental Techniques 


The experiments are performed using a Supersonic 
blow-down Wind Tunnel available at Birla Institute of 
Technology, Mesra, India. The test section was of the size 
of 50 x 100 mm giving a freestream Mach number of 2.0. 
Flow into the tunnel was established using a pressure 
regulated solenoid valve of 20 cm diameter. Settling 
chamber pressure of 3.08x10° N/m? was maintained 
which corresponds to freestream Reynolds number of 
3.6x10’ per meter. The pressure in the settling chamber 
was measured using a pressure transducer (Make Sensym, 
Model ASCXI50DN). Tests have been made on a cylin- 
drical protrusion with vortex generator placed ahead of 
the protrusions. The cylindrical protrusion had a diameter 
(D) of 24 mm and height (H.) of 3, 6 and 9mm. The ramp 
type vortex generator model has a height (Hyg) of 3mm, 
chord length of 12 mm, and ramp angle of 12 deg. The 
vortex generator was placed at a distance of 2.5 times the 
diameter (D) ahead of the protrusion. To capture the 
overall flow field, standard schlieren flow visualization 
technique was adopted. Oil flow visualization tests were 
conducted to obtain the surface flow patterns. Streak of 
flow on the surface plate was obtained by making use of 
a suitable mixture of Titanium dioxide, Oleic acid and 
lubricating oil. The oil flow patterns and schlieren flow 
photographs were captured using a digital camera (Model: 
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Sony DSLR A100K). Static pressures were measured at 
15 locations ahead of the protrusion by providing pres- 
sure ports of 0.8mm diameter on the surface plate, which 
were suitably connected to a 32 channel Electronic Pres- 
sure Scanner (Make: SCANCO, Model: ZOC 22B/32Px). 
Data was acquired using NI DAQ, LabVIEW software 
and a PC based Data Acquisition System. The details of 
the experimental set up are given in Ref [14]. The summary 
of the experimental test cases are tabulated in Table 1. 


Numerical Techniques 


Three dimensional structured grid was created using 
ICEM CFD. The dimensions of the grid were based on 
the diameter of the cylinder (D). The length, width and 
height of the domain are 35D, 12.5D and 12.5D respec- 
tively. The schematic of CFD domain along with boun- 
dary condition is shown in Fig.l. The close up view of 
grid over the model is shown in Fig.2. The origin was 
fixed at the center of the cylinder. Although the vortex 
generator geometry was simple, the computational grid, 
which fills the immediate surrounding volume, can be 


Table 1 Summary of experimental test cases 


H., mm H/D H./ ò 
Case 1 3 0.125 1 
Case 2 6 0.250 2 
Case 3 9 0.375 


Vortex generator : Hyg = 3 mm, 
c= 12 mm, Ap = 12° 


Case 4 
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Fig. 1 Schematic of CFD domain with boundary condition 


Fig. 2 Close-up view of grid around cylinder and cylinder 
with VG 
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challenging to build, especially for one using a conven- 
tional structured hexahedral grid generator. The triangu- 
lar shapes and sharp corners can lead to highly skewed 
cells adjacent to the device. In addition, the oblique in- 
duced shock can be difficult to resolve because it is 
usually not aligned with either family of grid lines. Many 
investigators use hybrid grids and/or over-set grids to 
overcome these difficulties [13]. 

In this study, Icem-CFD was used to produce multi 
block structured grid. The first cell size of the grid adja- 
cent to the bottom wall was Ay = 0.001 having wall y+ = 
0.26. The overall cell count in the flow domain for cy- 
linder was 100x75x65 and for cylinder with vortex ge- 
nerator was 100x75x85 approximately. The base of the 
computational domain was taken as no-slip wall boun- 
dary condition. Pressure inlet and pressure outlet boun- 
dary conditions were adopted for inflow and outflow 
boundaries. The top and side stream was taken as pres- 
sure far-field boundary condition. Walls of the cylinder 
and VG are considered as no-slip wall boundary condition. 

Compressible three dimensional numerical simula- 
tions were made to obtain the flow field around cylin- 
drical protrusions and with vortex generator at various 
Mach number using commercial CFD solver Ansys 
FLUENT. A density based coupled implicit formulation 
was used to solve the general equations of continuity, 
momentum and energy along with equation of state of 
gas. For the present computation k-œ turbulence models 
was used. The standard k-œ turbulence model uses two 
transport equations to solve turbulence parameters along 
with the general equation. Fluent uses finite volume me- 
thod ideally suited for both compressible and incompres- 
sible flows. The 3D multi block structured mesh was 
created around the cylindrical protrusion with VG. The 
Advection Upstream Splitting Method (AUSM) was used 
for solving conservation equations. This flux splitting 
method has an advantage over ROE flux splitting in 
terms of capturing of shock and discontinuities. The go- 
verning equations were solved using the second-order 
upwind scheme for convection terms including flow, 
turbulent kinetic energy and specific dissipation rate. It 
was able to capture most of the flow features around the 
protrusion. Ideal gas was used as the working fluid and 
its viscosity was calculated using Sutherland’s viscosity 
law. The courant number and Relaxation factor was in- 
itialized with 0.5 and 0.4 respectively for better conver- 
gence and stability at initial stage. Later they were in- 
creased gradually to a value of 2 and 0.8 for smooth 
convergence. For solution, convergence level of residuals 
was brought below 10-6 and simultaneously mass flow 
rate and drag coefficient value were also monitored. 


Results and Discussion 


The oil flow and schlieren flow visualization tech- 
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niques were used to study the flow features around the 
cylindrical protrusion at Mach number 2.0. Oil flow vi- 
sualization was done by using a suitable mixture of TiO2 
with oleic acid and SAE-40. Very fine drop of this care- 
fully prepared oil flow mixture was sprayed all over the 
surface plate around the protrusion and tests conducted 
for different height of protrusion with VG were forcibly 
drifted away along the free stream flow direction result- 
ing in a pattern of streak lines of the flow. The occurrence 
of shock waves ahead of the cylindrical protrusion was 
clear and distinctly observed. The comparison of experi- 
mental surface flow pattern and the numerical surface 
flow pattern obtained at a height of 0.01mm is shown in 
Fig.3. The region in the vicinity of the shock shows an oil 
flow pattern, which suggests an existence of a flow sepa- 
ration zone ahead of the protrusion for all three cases. 


Fig. 3 Comparison of surface flow pattern for case 1 


The reattachment takes place almost 1D distance 
downstream of protrusion. In cases with VG, the reat- 
tachment is observed further downstream at around 1.5D. 
On either sides of the model, a lateral separation of the 
flow is also witnessed. 

Schlieren flow visualization images for the cylindrical 
protrusion of case 3 with and without VG is shown in 
Fig.4. Schlieren images of cylindrical protrusion with 
single vortex generator clearly indicate that the adoption 
of vortex generator will improve the separated flow field 
ahead of the protrusion. The streamline pattern on the 
symmetry plane with presence of VG ahead of cylindric- 
al protrusion of case 2 is shown in Fig.5. The counter 
rotating vortices from the VG energize the boundary 
layer. The reattachment point moves further downstream. 
The comparison of pressure distribution ahead of the 
protrusion is shown in Fig.6. A good agreement between 
computed and measured pressures was obtained for all 
the cylinders of different heights. The increase in height 
of protrusion causes the upstream separation distance to 
increase. For the small height protrusion, even though the 
upstream separation zone is smaller when compared to 
the larger height protrusion, the flow passing over the 
protrusion affects the downstream flow field. 
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(a) (b) 


Fig. 4 Schlieren flow visualization (a) case 3 and (b) case 3 
with VG 
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Fig. 5 Streamline pattern around symmetry plane (a) case 2 
and (b) case 2 with VG 
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Fig. 6 Measured and computed Cp for cylindrical protrusion 
with different heights 


The difference in peak pressure observed from Fig.6 
also indicates the decrease in pressure with respect to the 
protrusion height. The comparison of case 3 protrusion 
with and without VG is shown in Fig.7. Case 3 protru- 
sion with VG shows a good reduction in peak pressure by 
about 36%, which can be advantageous. 

Figure 8 represents the pressure distribution on cylin- 
drical protrusions of various heights along the centerline 
(Z/D = 0) with and without a vortex generator. The VG 
also helps in reduction of pressure acting on the cylin- 
drical protrusion. From Fig.8, it is observed that the 
pressure distribution on cylindrical protrusion with VG 
follows similar trend. 
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Fig. 7 Computed Cp for case 3 protrusion with VG 
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(b) at X/D = 0.5 and Z/D = 0 


Fig. 8 Stagnation line pressure distribution for different height 
cylinder with VG 


The drag acting on the cylinder is estimated numeri- 
cally using Fluent. It is observed that the aerodynamic 
drag acting on the cylinder is reduced by adopting VG 
ahead of the cylinder. The comparison of drag coefficient 
of different height cylinders with VG is shown in Fig.9. It 
is seen that around 30% reduction of drag is achieved 
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through adopting VG ahead of protrusion of case 3 and a 
reduction of around 15% and 8% is observed for case 2 
and case | respectively. This is because of the counter 
rotating vortices which energize the boundary layer se- 
paration and also increase the pressure distribution on the 
rear surface of the cylindrical protrusion. This in turn 
leads to an overall increase in downward force. 
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Fig. 9 Computed drag coefficient for different height protru- 
sion with VG 


Conclusions 


The present study was carried out to obtain a clear 
view on the flow field around a small cylindrical protru- 
sion which is of the order of the local boundary layer 
thickness. The effect of the height of protrusion on the 
flow field were investigated experimentally and numeri- 
cally for a supersonic flow at Ma =2.0. The presence of a 
complex flow field around the protrusion due to shock 
wave boundary layer interaction was observed from flow 
visualization techniques. A flow separation is observed 
ahead of the cylindrical protrusion due to the protrusion 
bow shock and surface boundary layer. This flow separa- 
tion produces a separation shock which coincides with 
the protrusion bow shock and creates a lambda shock 
structure. The primary separation line that moves up- 
stream with increasing height of protrusion is also seen 
from the plots. The reduction in peak pressure ahead of 
the cylindrical protrusion is observed with the presence 
of VG. The VG helps in improving the separation region 
ahead of protrusion by energizing the boundary layer 
downstream. The reduction in pressure drag on cylin- 
drical protrusion was achieved by adoption of vortex 
generator ahead of cylindrical protrusion, which im- 
proves the overall aerodynamic performance. 
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